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— ». I We study the vibration modes of the Skyrme model within the rational map 

^^ ' ansatz. We show that the vibrations of the radial profiles and the rational maps 

are decoupled and we consider explicitly the case B = 1, B = 2 and B = 4. We 

then compare our results with the vibration modes obtained numerically by Barnes 

QJ \ et al. and show that qualitatively the rational map reproduces the vibration modes 
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obtained numerically but that the vibration frequencies of these modes do not match 
very well. 
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Proposed by Skyrme as a fundamental theory of strong interactions, the Skyrme[T] model 
was later shown by Witten[2j to be a low energy limit of QCD in the limit of large 
r> I colour number. In that context, the classical solutions of the Skyrme model correspond 

c^ ' to bound states of QCD. The simplest solution, with baryon number 1, can be computed 

analytically up to solving an ordinary differential equation numerically. For larger baryon 
numbers, one must compute the solutions by solving the full classical equation of the 
model numerically [3]. 

Recently, Houghton et al.[l] showed that the solution of the Skyrme model can be 
well approximated by the so called Rational Map ansatz. In this ansatz solutions are 
approximated by a radial profile function and a rational map ansatz which only depends 
on the polar angles variables. One then determines the radial profile by solving an ordinary 
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The configurations obtained by the ansatz fit the numerical solutions very well. 

Once a classical solution has been obtained, one must still quantise some of the re- 
maining degree of freedom. One way to do this is to compute the vibration modes of these 
solutions. This was done for the numerical solutions of the model by Barnes et al. ^ ^ 
for baryon numbers B = 2 and 4. In this paper we compute the vibrational mode of the 
rational map configuration and compare them to those obtained numerically. 

The Skyrme Model [l] is defined by the following Lagrangian 

+^^Tr{UP-l)^d^x, (1) 

where U = U{x,t) is an SU{2) chiral field, F^^ is the pion decay constant, 771^^ is the 
pion mass and e a parameter of the model which is determined by fitting the classical 
solutions to experimental data. Notice that the last term in ([T]) is the so called mass term 
where we have used the generalised mass term proposed by Kopeliovich et al. [7| where 
the parameter p is a positive integer. 

Rather than using ([T]), it is convenient to rescale the space-time coordinates and the 
mass parameter, x^ = 2x^/F.„e; m = 171^,2/ F.„e and use the dimensionless Lagrangian 

e 
' - Lp 



Stt^F, 



^2^2 J I - iTriR.nn + ^Tr{[R„ i?.][i?^ i?l) + ^Tr{U^ - 1)} rf^x , (2) 

where i?^ = {d^U)Ul 

To approximate the solution by the rational map ansatz, we first introduce the complex 
coordinate ^ = tan(|)e*'''' where 9 and 4> are the polar angles. Then the rational map 
ansatz is given by 

U = ^2if{r)n^^^ya ^ (3) 

where a are the Pauli matrices and 

nRio = 1^ 1^(^)1 2 {2^{R{0),2Q{R{0)A-\Rm') ■ (4) 

The degree of the rational map R{^) corresponds to the baryon number of the configura- 
tion (see [1]). 

To approximate the classical solution of a given baryon number B, one would thus 
take 
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Knowing the value of I, one then uses the Euler-Lagrange equation to derive the equation 
that the profile /(r) must solve. In [^, it was shown that the case S = 1 is nothing but 
the hedgehog solution computed by Skyrme. For B = 2, the configuration is axially 
symmetric while for i? = 4 it has the symmetry of a cube. 



2 Vibration Modes 



To study the vibration modes of the rational map ansatz configurations minimising ([2]), 
we add a time dependant perturbation to the rational map and the profile function around 
their minimising values. We then insert the perturbed ansatz into the Lagrangian ([2]) and 
compute the Euler-Lagrange equation for the perturbation, keeping only the linear terms. 

Denoting fo{r) the minimising profile function for the static solution, we take 

f{r,t)=fo{r)+gir,t) (7) 

where g is assumed to be a small fluctuation around /o satisfying the boundary condition 
giO,t)=g{oo,t) = 0. 

To perturb a rational map, we must perturb all its coefficients, even the one that are 
null, by adding a small time dependant perturbation. 



R{^,t)=P{U)/Q{U) 



where 



P{^, t) = Ei=o^(a. + 6a,)e = PoiO + SP{^, t) 

Q(e, t) = EiZoibj + 6bj)e = QoiO + m^, t) (9) 

Notice that, as i? is a ratio of P and Q, the coefficients of the rational map are determined 
up to an overall constant. If as is non zero, we can divide both P and Q by aB{l+SaB/aB) 
to linear order and Sqb then can be incorporated into the other 6ai and Sbi] (if as is null, 
one divides P and Q by 6^(1 + Sbs/bB))- The rational map perturbations are thus 
described by 4(5 + 1) — 2 parameters. 

When inserting (^ into ([2]), the perturbed rational map only occurs in the integral 
([6]) and in the two expressions 

^ ^ j_ f \R? 2zd^d^ 



An J (l + |i?|2)2(l + |^|2)2 ' 

^ 1 r (l + I^P)^ cii^2 2zd^d^ 

^ Air J (l + |i?|T d^ (l + I^P)' ■ 

3 



tto = i^or+iQor 

ai = Poi6P) + Poi6P) + Qoi5Q) + Qo{5Q) 

«2 = i6P)i5P) + {6Q){6Q) 

/3o = Po,(Qo — PoQo,^ 

A = Po,dSQ)-Po{6Q)^-Qo,dSP) + QoiSP)^ 

P2 = i5P)^iSQ) - iSP)iSQ)^ 

70 = m' 

71 = 2\Po\'Wi+0o/3i) 

Ai = QoiSP) - Po{6Q) 



;ii) 



121 



we have 



I 



1 

47r 



;i + ieiY 



On 



+ (- 



7o+ 7i 



470^) 



3^ 



1 

47r 
1 



+ 107o4 - 470 471— + 72 

^ Oq ao ao 

Alp 2id^d^ 



2ididi, 



al 



;i + iei 



2^2 



«0 



2i di di , 



Notice that the integrals (EI) and ( IT3l) are at most quadratic in the parameters 
To rewrite these integrals in matrix form, we define 
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and rewrite (IT3l) as 



I 

y 



rT-- 






(15) 



where Xq is the value of I for the unperturbed rational map as given in [Jj. Notice that 
there is no linear term in V . This is because the unperturbed rational map minimises (EI) 
and (|TU|). 



Inserting the perturbed ansatz ([7]) and (jH]) into ([2]), we get 



L = Lo + L2 (16) 



where 



Lo = ^|{-25sinVo-(25smVo + 0/o'-^sinVo 
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Notice that the perturbation of the radial profile and the rational map are completely 
decoupled. We can thus study the radial and angular vibration modes separately. To do 
this we must compute the Euler-Lagrange equations for g and for V to linear order and 
solve the resulting eigen value equations. 



3 Rational Map Vibrations 



The equation for the rational map vibrations V is straightforward to derive from flTUj) and 
can be written as 

A,,V, = -D,,V, , (18) 



where 



Aij — X2ij Ti + y2ij T2 

A, = X2.,r3 (19) 



and 



T2 = j (sin"^ /o) dr 

Ta = J{^,sm'fo)dr (20) 

are numerical parameters which can be evaluated numerically. We provide their value for 
the case m.,^ = and m,r = 0.526 and p = 1 in Tables [T] and [21 
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of the entries can be shown to vanish using parity symmetries; we have evaluated the 
others using both Maple and Mathematica. 



Then we inserted V = Vosin(ct;t) into (fTSl) and used Maple and Mathematica to solve 
the resulting eigen value problem. We would like to point out that the eigen vectors, i.e. 
the vibrations, do not depend on the actual values on Fi, r2 and Fa (or p and m), while 
the eigen values, i.e. the vibrations frequencies do. The results are described separately 
for the cases B = 1, B = 2 and B = 4. We will then compare our results with those of 
Barnes et al. [5] [6] obtained, in our conventions, for the value m.,r = 0.526. 



3.1 B=l 



For the hedgehog solution the unperturbed ration map is simply R = ^ and there are 6 
perturbed rational maps: 



Rbtr,x{C) = ]Tj:d|; RbtT,y{C) = l+idV ^btr,z{0 = T 
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with 



d = ecos{ujt) 
where e and u are respectively the (small) amplitude and the frequency of vibration. 



(21) 
(22) 



Rrot,x, Rrot,y and Rrot,z are zero modes and correspond respectively to a rotation around 
the X, y and z axis. Rbtr,x, Rbtr,y and Rbtr,z have the same eigen value uur = 0.7575 for 
171^^ = and uJbtr = 0.9239 for m^r = 0.526. They form a 3 dimensional eigen space as they 
are conjugated to each other through a 90° rotation. This eigen subspace corresponds 
to the broken translational invariance of the solutions: while the Skyrme Lagrangian 
([2]) is invariant under translation, the rational map ansatz breaks that symmetry by 
pinning the center of the solution. If one performs the translation x -^ x + Xq, into the 
hedgehog ansatz, the ansatz is broken, as the radial profile becomes a function of the 
polar angles and the rational map becomes a function of r. If xq is small, one can expand 
the translated expression, keeping only the linear terms in xq. The rational map then 
becomes -R(0 = i+esJr ■ ^^ ^^^ energy density of the Skyrmion is concentrated around 
a shell of radius tq, we can take R{$,) = i,^cj° as an approximation of the translated 
rational map, recovering Rbtr,x{0- -^^ ^^ had to make several approximations to derive 
that expression of the perturbed rational map, instead of being a zero mode, it has a non 
zero vibration frequency that does not have any physical meaning. 



B= 1 
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4.28869 

0.872418 

0.370159 



7.58651 
0.96429 
0.16419 



12.2868 
1.02577 
0.0895335 



Table 1: Fi, F2, F3, for m^r = 



does not have any genuine non-zero vibration mode. 



3.2 B=2 



The unperturbed rational map for the case B - 
perturbed rational maps. 

First we have the 3 rotational zero modes 



2 is given by i? = i^^ and there are 10 
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Then there are 2 iso-rotation zero modes 
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(24) 



There are only 2 such modes because the iso-rotation around the z axis coincides with 
the proper rotation around the z axis. 



The lowest vibration modes are 
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(25) 



but they correspond to the broken translation modes. Their eigen values uur are respec- 
tively 0.69, 0.69 and 0.4766 for m^ = and 0.847 0.847 and 0.584 for m^ = 0.526, but 
they have no physical meaning. 



We are thus left with a 2 parameter class of genuine vibration modes: 
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These two perturbed rational maps correspond to lateral squeezing and stretching of the 
torus alternating between elongations along the x and y axis as shown in Figure [H This 
mode is sometimes referred to as the scattering mode as it corresponds to two Skyrmions 
colliding with each other. Rs2 corresponds to the same deformation but rotated by 45°. 
By taking a linear combination of Rsi and Rs2 the torus can be made to wobble along 
any axis in the x — y plane. The eigen value for these 2 modes is Ug = 0.9909 for m^ = 
and cUo = 1.2255 for m^ = 0.526. 
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3.04713 
0.75611 
0.405799 



5.40344 
0.82419 
0.185219 



8.94968 

0.885387 

0.103598 



Table 2: Fi, F2, F3, for mj, = 0.526 and p = 1 





Figure 1: B = 2, Rsi{0- (Left) e > 0; (Right) e < 0. 



3.3 B=4 



The unperturbed rational map for the case 5 = 4 is given hj R = Po/Qo where Pq = 
^^ + 2-\/3 i ^^ + 1 and Qo = ^"^ ~ 2\/3 i ^^ + 1. There are 18 perturbed rational maps, 
including 6 zero modes: 
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Pq + (22 9, - 2^3 9,)e + (2z 9, - 2^3 g^)^ 
Qo + (2z ^:. + 2^3 ^,)e + (2z ^, + 2^3 9,)^ ' 
Pq + (2v^ i g^ - 29y)e + (-2v^ i g^ + 29y)^ 
Qo + {-2V3 1 9y - 29y)i^ + (2^3 1 9y + 29y)i ' 
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which correspond to the 3 rotations Rrot,x, Rrot,y and R 

-^iso.y anQ -thiso^Z' 

The lowest vibration modes are the broken translation 
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The first genuine vibration mode is given by 
Rtb,y{0 

Rtb,ziO 
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and its eigen value is ojtb = 0.6093 for m-,, = and utb = 0.7508 for m^r = 0.526. To picture 
it one, fias to think of the cube as two tori stacked on top of each other. Then each torus 
oscillates like the scattering mode of the B = 2 solution but in phase opposition. This 
results in a cube which is somewhat squeezed into a toroidal shape as shown in figure 
El This can also be thought of as the scattering of four Skyrmions in a tetrahedral 
configuration. 





Figure 2: S = 4, RthAO- (Left) e > 0; (Right) e < 0. 

The last triplet mode corresponds to the pinching and stretching of two opposite edges 
of the cube so that two of the opposite faces of the cube are deformed into a rhombus. 
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(30) 



The energy density plots are presented in figure The eigen value of the rhombus mode 
is oJrhomb = 0.7395 for m^r = and oJrhomb = 0.9113 for ttIt, = 0.526. 
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Figure 3: B = A, Rrhon^bAO- (Left) e > 0; (Right) e < 0. 



The doublet mode 



Rscatl (0 
Rscat2\C) 



Qo + ie^^ + ie' 
Po + iee 



(31) 



corresponds to a deformation where the cube is ahernatively stretched and then flattened 
along one of the symmetry axis going through the center of the cube's faces. A combi- 
nation of a stretch along two of the (perpendicular) axis is equivalent to a stretch along 
the third axis. This corresponds to the scattering mode of two tori along their axis of 
symmetry as shown in figure HI The eigen value for this class of modes is Ustr = 0.818 for 
itIt, = and ujstr = 1.0064 for m^ = 0.526. 





Figure 4: S = 4, Rscat2iO- (Left) e > 0; (Right) e < 0. 
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R. 



tet 



Qo(l-e) 



(32) 



corresponds to a tetrahedral deformation of the cube as illustrated in Figure O Its eigen 
value is Utet = 1-1355 for m^ = and Ufet = 1.4026 for m^ = 0.526. 





Figure 5: 5 = 4, RtetiO- (Left) e > 0; (Right) e < 0. 



4 Radial Vibrations 



The Euler-Lagrange equation for the radial perturbation g can be easily derived from 
(HE]): 

{■g - g") (2B sin^ f, + r^) - (2r + 2Bf, sin(2/o)) g' - (-2B cos(2/o) + 2i3/^2 ^^^(3/0) 



+2S/^' sin(2/o) - ^ [6(sin2 /„) (cos^ /o) - 2 sin^ h] - m'r^ cos(p/o)) ^ = 0, (33) 



-1 /o) (cos /o) - 2sm /oj - m r 

where g' = -^ and g = -^. Taking a perturbation of the form 

9ir,t) = goir) sm{ujt) (34) 

and inserting it into fl33|) leads to a Sturm-Liouville equation where a;^ is the eigen value. 

Notice that when uj > m, any perturbation is radiated away and there are no genuine 
vibration modes. When the pion mass 777,7^ is too small, there are thus no genuine radial 
vibration modes. On the other hand, the Skyrmion still has so-called pseudo-vibration 
modes, as studied by Bizon et al. p], where the excitation energy is radiated away relatively 
slowly. In a quantum theory these modes would correspond to resonances. 



We have solved (IHHjl numerically for several baryon numbers and several values of the 
pion mass and the parameter p. Our results are summarised in figure E] where we present 
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Figure 6: Radial Vibration frequency for i3=l,i3=2, i3=4. 
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±J \ p 










1 


1.60 


1.19 


0.96 


0.85 


2 


1.10 


0.82 


0.71 


0.65 


4 


0.77 


0.59 


0.53 


0.51 


6 


0.60 


0.47 


0.44 


0.43 


7 


0.56 


0.44 


0.42 


0.41 


9 


0.47 


0.38 


0.36 


0.36 


17 


0.33 


0.28 


0.27 


0.27 



Table 3: Critical Mass as a function of B and p. 

the vibration frequency cu as a function of the pion mass for different baryon number 
values. 

On table [S] we present the critical value of the pion mass below which there is no 
genuine vibration mode. Notice that the case 5 = 1 is the exact hedgehog solution and 
corresponds to a single nucleon. The critical value of the mass for that case, can thus 
be interpreted as an upper bound for the pion mass when it is taken as an adjustable 
parameters [9] [10], as the proton has no stable excited state. 



5 Comparison 



Several years ago, Barnes et al[5][6] computed the vibration modes of the B = 1 and 
i? = 4 Skyrmion solutions numerically. In their work, they only considered the standard 
mass term, p = 1 and the pion mass that they used, in our parametrisation is given by 
m^r = 0.526. By comparing their results with ours, we will be able to asses the quality of 
the rational map for the estimation of the vibration modes. 

The method used by Barnes et al. consists in solving numerically small fluctuations of 
the Skyrme field around a static solution obtained numerically too. To solve the Skyrme 
equation numerically, one has to invert a matrix that is field-dependant. To improve the 
efficiency of their code, and motivated by the fact that they only studied small vibrations, 
Barnes et al. set that matrix to the one obtained from the stationary field. This effectively 
pinned the position and iso-orientation of the Skyrmion, thus breaking the translation and 
iso-rotation zero modes. On the other hand, the rotations symmetry was preserved, apart 
from the symmetry breaking introduced by the square lattice and the periodic boundary 
conditions. 

We compare our results with those of [3] and [U] in Tables H] and One can identify 
the various vibration modes by comparing the dimension of the subspace they span, i. e. 
the dimension of the representation they belong to, as well as their symmetries and the 
deformation induced in energy density plots. Figures 2 in [5j match exactly those we have 
obtained for the vibrational modes of the rational map of -B = 4. Figure 3 in [6j, once 
translated to what is plotted in figure [T], also matches our observation for B = 2. 

For B = 2, we had to read the values of u from figure 1 in [6j. There is only one 
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rational map for that mode is much higher than the one obtained numerically, but the 
frequency obtained from Barnes et al. indicates that it is a genuine vibration mode. 

The mode described as a dipole breathing motion in ^, turns out to be a broken 
translation. Notice also that the breathing mode is a pseudo vibration mode for both the 
numerical and the rational map results. The radiative modes observed in [B] are the results 
of the numerical method used and are excluded from the rational map ansatz, except from 
the mode with symmetry Aiu which happens to be a broken translation mode. 

For B = A, the modes that we obtained match those obtained in [Sj. The mode with 
^num = 0.655 in [5] turns out to be a broken translation mode rather than a proper 
vibration mode. This is supported by comparing the energy densities for that mode, not 
presented in this paper, and the one given in [5], as well as the symmetry of that mode. 
Apart from that, the rational map ansatz predicts 5 non-null vibration modes out of 6 
(the diagonal mode, with Unum = 0.738, is missing). Symmetry considerations suggest 
that the radiative mode with Unum = 0.587 might be a pseudo-iso-rotational mode, but it 
is difficult to be conclusive as one would need to know more about this numerical mode to 
make a definite statement. If the vibration modes obtained from the rational map ansatz 
match the one obtained numerically rather well, the predictions for the frequencies of 
these modes are rather poor and are not even in the correct order. Overall, the rational 
map ansatz appears to be stiffer than real Skyrmion solutions. 



6 Conclusion 



The rational map ansatz has been shown to be a good approximation to the shell-like 
solutions of the Skyrme model, i.e. when the pion mass is null or relatively small, or 
when the baryon number is small. In this paper we have shown that the rational map 
ansatz can also be used to study the vibration modes of the Skyrmion solutions and that 
qualitatively it produces the correct vibration modes. For i? = 1, for which the ansatz 
is actually the exact solution, and B = 2, the vibration modes are all predicted, but for 
B = 4 they are all predicted except for one, the one with the largest vibration frequencies. 

On the other hand, the vibration frequencies obtained from the rational map ansatz 
are quite poor when compared to the frequencies obtained numerically in [H] and [H] . The 
relative order is not even correct. We interpret this result by saying that the rational map 
ansatz is too stiff and that the configuration is restrained to vibrate in a subspace that is 
too narrow. 

Our study has also shown that some of the modes observed by Barnes et al. were 
mistakenly taken as genuine vibration modes when they are actually broken vibration 
modes. 
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or 



Numerical [6] 



Ansatz 



UJ 


Degeneracy 


Description 


Symmetry 


LU 


Degeneracy 


Description 


Symmetry 


0.03 


2 


broken rotational modes 
(around x and y axes) 


E\g 





3 


X, y rotational modes {Rrot,x,Rrot,y) 
and z (iso)rotational mode {Rrot.z) 


Elg + A2g 


0.31 


2 


2 Skyrmions scattering mode 


E2g 


1.225536 


2 


2 Skyrmions scattering mode {Rsi,Rs2) 


E2g 


0.49 


1 


radiative mode 


Alu 


















0.583500 


1 


broken z translation {Rbtr,z) 


Alu 


0.52 


2 


radiative mode 


E2g 










0.75 


1 


breathing mode 


Alg 


0.76 


1 


breathing mode 


{Alg) 


0.84 


2 


a dipole 'breathing' motion^ 


Elu 


0.846727 


2 


broken x, y translation {Rbtr,x, Rbtr,y) 


Elu 



Table 4: Comparison of vibrational frequency between the numerical results and the rational map ansatz, B = 2. 

' Its description of the energy density matches that of the broken x, y translational zero mode obtained from the rational map ansatz; however, this 
representation, Eiu, coincides with the doublet of x, y translational zero mode of the 2-monopole toroidal BPS solution. 



Oi 



0.070 
0.367 
0.405 
0.419 
0.513 

0.545 

0.587 
0.605 
0.655 
0.738 
0.908 



Numerical fsl 



Ansatz 



Degeneracy Description 



broken rotation 
deuteron scattering mode 
tetrahedral mode 
rhombus mode 
4 Skyrmions mode 

radiative mode 



nimctry 


U! 







^13 





E, 


1.006370 


A2„ 


1.402608 


F2, 


0.911260 


F2„ 


0.750841 



radiative mode 




A2u 




breathing mode 




Alg 


0.64 


broken translation 




Fi^ 


0.562538 


diagonal mode 




F2, 




lowest nonzero radiative 


mode 







Degeneracy Description 



iso-rotation {B.iso.x , Riso.y) 

rotation {Rrot,j2, Rrot,y, Rrot,z) 

deuteron scattering mode {Rscati-, Rscat^) 
tetrahedral mode (Rtet) 

rhombus mode {Rrhomb,x : Rrhomb.y : Rrhon 

4 Skyrmions mode {Rtb,a: , Rtb^y , Rtb,z ) 



iso-rotation {Riso,: 



>,.) 



breathing mode 
broken translation {R^ 



5 Rbtr.y, Rbtr,z) 



Table 5: Comparison of vibrational frequency between the numerical result and the rational map ansatz, S = 4. 

' Its description of the energy density matches that of the broken x, y translational zero mode obtained from the rational map ansatz; however, this representation, 
Eiu, coincides with the doublet of x, y translational zero mode of the 2-monopole toroidal BPS solution. 
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